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ABSTRACT 

Several young supernova remnants exhibit thin X-ray bright rims of synchrotron radiation at their 
forward shocks. Thin rims require strong magnetic field amplification beyond simple shock compres- 
sion if rim widths are only limited by electron energy losses. But, magnetic field damping behind 
the shock could produce similarly thin rims with less extreme field amplification. Variation of rim 
width with energy may thus discriminate between competing influences on rim widths. We measured 
rim widths around Tycho’s supernova remnant in 5 energy bands using an archival 750 ks Chandra 
observation. Rims narrow with increasing energy and are well described by either loss-limited or 
damped scenarios, so X-ray rim width-energy dependence does not uniquely specify a model. But, 
radio counterparts to thin rims are not loss-limited and better reflect magnetic field structure. Joint 
radio and X-ray modeling favors magnetic damping in Tycho’s SNR with damping lengths ~1 5% of 
remnant radius and magnetic field strengths ~50-400 fiG assuming Bohm diffusion. X-ray rim widths 
are ~1% of remnant radius, somewhat smaller than inferred damping lengths. Electron energy losses 
are important in all models of X-ray rims, suggesting that the distinction between loss-limited and 
damped models is blurred in soft X-rays. All loss-limited and damping models require magnetic fields 
> 20 /xG, affirming the necessity of magnetic field amplification beyond simple compression. 

Keywords: acceleration of particles - ISM: individual objects (Tycho’s SNR) - ISM: magnetic fields 
— ISM: supernova remnants — shock waves — X-rays: ISM 


1. INTRODUCTION 

Electrons accelerated in the forward shocks of young 
supernova remnants (SNRs) emit synchrotron radiation 
strongly in the shock’s immediate wake at radio wave- 
lengths and sometimes in X-rays. In a few cases, they 
quickly turn off downstream, producing a shell-like mor- 
phology of bright X-ra y and radio rims /filaments due to 
line-o f-sight projection (Bamba et al. 2003 Reynoso et al. 
1997 1 . Strong and time-var i able synchrotron radiation 


e.g., Uchiyama et al. 2007 Patnaude & Fesen 2007), 


in conjunction with mu ltiwavclcngth spectra l modeling 


( Aharonian et al. 2004 Acero et al. 2010 Ackermann 
et ai. 2013), suggests that electrons are accelerated to 
TeV energies in young SNRs. Although synchrotron 
emission due to accelerated electrons does not require ac- 
celeration of an unseen hadronic component, the prevail- 
ing theory of diffusive shock acceleration (DSA) should 
operate on both positive ions and electrons. Efficient 
hadron acceleration in supernova remnant shocks is a 
prime candidate source for galactic cosmic rays up to 
the cosmic ray spectrum’s “knee” at around 3 PeV (Vink 
2012). But many fundamental questions about shock ac- 
celeration remain unanswered. Under what conditions 
do shocks accelerate particles efficientl y? How aremag- 
netic fields amplified in such shocks? Reynolds ( 2008 ) re- 
views relevant observations and open questions to date. 
These questions are relevant to man y astrop hysi cal set - 


tings, such as Ear th’s bow shock ( Ellison ct al. 1990) 


starburst galaxies (Heckma n et a l. 
galactic nuclei (Ichen et ai.||20I4|), 


1990|, jets of 


active 


galaxy clusters (van 
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Weeren et al.||2010|), and cosmological shocks (|Ryu et al. 

DDSf 

Spectral and spatial measurements of synchrotron rims 
can constrain downstream magnetic field strength and 
structure. If rim widths are set by electron energy losses, 
post-shock magnetic fields must be amplified to ~10 2 rG 
to account for the thinness of observed rims ( Vink fe 


Laming||2003| [Bamba et al.||2003[ | 2005 [ [Voik et aT" 

the 


Parizot et ai.| 2006 ) . In these models, the magnetic fief' 
is assumed advected downstream and nearly constant 
over rim widths. Alternately, the magnetic field strength 
may be damped downstream of the shock and prevent 
electrons from radiating efficiently, so that thin rims re- 
flect magnetic field variation rather than efficie nt par- 
ticle acceleration and synchrotron cooling (Pohl et al. 


2005 1 . Damping, in particular, may permit less extreme 


magnetic field amplification. We refer to these as “loss- 
limited” and “damped” models for rim widths. We shall 
see that models can range continuously between these 
two cases, and that the distinction between the two can 
vary with observing frequency. 

The possi bility o f dam ping in S NR sho cks has not been 
fully tested. Marcowith & Casse (2010) compared phys- 
ically motivated magnetic damping models to X-ray rim 
widths and synchrotron spectrum cut-offs and thus sug- 
gested that only young SNRs (age < 500 yr) can exhibit 
magnetic damping in con junc tion w ith efficient particle 
acceleration. Rettig & Pohl (2012) gave model predic- 
tions for several historical SNRs and proposed discrim- 
ination based on filament spectra - the expectation is 
that damped spectra are softer, loss-limited harder. 

Hydrodynamic models can reproduce X-ray rim pro- 
files reasonably well with both loss-limited and damped 
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magnetic field m odels ( Cassam-C hcnai ct al.||2007 Mor- 


lino & Caprioli 2012 Slane et al. 2014). However, loss- 
limited models generally cannot reproduce thin radio 
rims. Radio-emitting GeV electrons do not lose substan- 
tial energy via radiation, so modeled intensities rise to a 
broad maximum toward the remnant interior, then grad- 
ually drop due t o spher icit y eff ects as the density drops 
in the interior. Reynolds ( 1988 ) empirically modeled ra- 
dio rims in the renmantT of SN 1006 and concluded that 
such thin rims required sharp gradients in electron energy 
density (from, e.g., time-v aria bl e partic l e accele r atio n) 
or magnetic field strength. Cassam-Chenai et al. (2007) 
used a 1-D hydrodynamic model with nonlinear L)SA to 
jointly model radio and X-ray rims in Tycho’s supernova 
remnant; neither loss- limited nor damped models could 
match radio profiles and radio/X-ray intensities simul- 
taneously. Moreover, gradual X-ray spectral variation 
observed downstream of the shock front was poorly re- 
produced with both models (Cassam-Che nai et al.|2 007). 
But damped models were able generate limb- brightened 
radio rims at the forward shock, even if morph ology was 
not en t irely c onsistent with observation. |Cassam-Chenai| 


et al. (2007) thus suggested that some combination of 
amplification and magnetic field variation might explain 
radio morphology. 


Recently, Ressler et al. (2014) (hereafter, R14h sought 
to discriminate between damped and loss-limited rims by 
measuring rim width-energy dependence in X-ray ener- 
gies in SN 1006. In the simplest models, rim widths are 
expected to be roughly energy-independent if rims are 
damped, whereas widths should narrow with increasing 
energy if rims are energy loss-limited. R14 included a 
variety of effects which can blur this distinction. To fur- 
ther test these models, we follow |R14| by measuring X- 
ray rim widths at multiple energies in Tycho’s supernova 
remnant (hereafter, Tycho). Tycho exhibits an exten- 
sive shell of synchrotron-dominated thin rims around its 
periphery (Figure [l]); the rims show very little thermal 
emiss ion, consistent with e xpansion into a low density 
ISM (|Williams et al.||2013h . A deep 750 ks exposure of 
the entire remnant from 2009 allows fine sampling of the 
remnant rims. We build upon previous estimates of mag- 
netic field strength and particle diffusion in Tycho that 
draw from multiwavelength obse rvations and various as- 
sumptions on CR acceleration ( Vojhetal. 20021 20051 
Parizot et al.||2006l |Morlino fe (JaprioTT 2012| Rettig fe 


Fold 20 1 2) e.gij 


Uur procedure cl osely follows that of |R14| We first 
review the model of |R14| used to model rim profiles and 
widths, then describe the procedure for selecting, mea- 
suring, and fitting rim widths to model width-energy de- 
pendence. We explore degeneracies in model fitting and 
consider radio rim morphology as an additional discrim- 
inant between models. 


2. NONTHERMAL RIM MODELING 
2.1. Particle transport 

The energy and space distribution of electrons at a 
supernova remnant’s forward shock controls the syn- 
chrotron rims we see in X-ray and radio. We assume 
that diffusive shock acceleration generates a power law 
distribution of electrons with an exponential cut-off at 
the forward shock and model 1-D steady-state plane ad- 
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Figure 1. RGB image of Tycho with region selections overlaid. 
Image bands are 0.7-1 keV (red), 1-2 keV (green) and 2-7 keV 
(blue). Bold region labels (1, 16) indicate region selections shown 
in Figures [3] [T] Filament 1: Regions 1-3, filament 2: regions 4-10, 
filament 3: regions 11-13, filament 4: regions 14-17, filament 5, 
regions 18-20. 


vection and diffusion of the electron distribution f(E , x), 
where E is electron energy and x is distance downstream 
of the forward shock: 

^£-1 ( D S)“i VE = XoE-e-^-m, 

111 


Following 

Berezhko & Volk 

2004); Cassam-Chenai et al. 

(2007 

|;M 

iorlino et al. (2010 

; Rettig &Fohl (20l2f. The 


ing downstream of the shock; D is the diffusion coeffi- 
cient, Vd is fluid velocity downstream of the shock, and 
the constant b = 4e 4 /9 m 4 c 7 = 1.57 x 10 -3 in appropri- 
ate CGS units arises from synchrotron power loss (i.e. , 
dE/dt = —bB 2 E 2 , averaged over pitch angles). The ini- 
tial electron distribution is specified by an arbitrary nor- 
malization Kq, DSA cut-off energy E cut (givenjn Sec- 
tion | 2.3 ), and spect ral index s = 2a + 1. Zirakashvili 


& Aiiaronian (2007) derive an electron energy spectrum 
with super-exponential cut-off e~^ p ^ Pcut ' > , but we use a 
simple exp onential cut-off for simplicity and consistency 
with R14| We have not yet specifie d th e cause of the 
injected - spectrum cut-off (see Section 2.4), but the func- 


tional form of a power law with exponential cut-off is a 
good approximation to predictions for DSA spectra lim- 
ited by synchrotron losses, remnant age, or p article es- 


cape ( Webb et al.||1984 Reynolds 1998 2008). All con- 


stants aiid - equatIons are given in CCS (Gaussian) units. 
Our presentation is somewhat abbreviated, but a fuller 
exposition and literature review are given by |R14| 

For Tycho , we a dopted radio spectral index a = 0.58 
(Sun et al. 2011) and hence electron spectral index 
s = 2a + 1 = 2.16. We assumed a remnant distance 


3 kpc (but cf. Hayato et al. 2010), which gives a shock 
radius of 1.08 x 10 l a cm from the observed angular radius 
240" (jGreen 2014) and further sets shock velocities for 
each rim profile we consider. Tycho’s forward shock ve- 


locity varies with azimuth by up to a factor of 2 (Katsuda 
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et ab 2008); we linearly interpolated velocities reported 
by Williams et ai. (2013) (rescaled to 3 kpc) to estimate 
individual shock velocities for each region. We assume a 
compression ratio of 4 as for a strong shock (unmodified 
by cosmic-ray pressure), and take downstream velocities 
Vd to be one-fourth the interpolated shock velocities. 

We assume isotropic diffusion and only consider par- 
ticle transport downstream of the forward shock. The 
velocity is assumed constant, as is the magnetic field 
for loss-limited model rims, in contrast to the expected 
Sedov-Taylor similarity solution for velocity and density 
in an adiabatic blastwave. Our assumptions of constant 
velocity, magnetic field, and plane flow should be reason- 
able as we generally consider synchrotron emission within 
10 % of the shock radius, r s , from the forward shock, 
though a few models are followed far enough inward that 
this approximation begins to break down. We consider 
profile emission strictly upstream of both th e conta ct dis- 
continuity and reverse shock (cf. Warren et al. 2005), 
and avoid regions where the contact discontinuity over- 
runs the forward shock. Our modeling neglects flow and 
electron spectrum modification due to the nearby con- 
tact discontinuity and particle acceleration at the reverse 
shock. These effects could matter (and we do not quan- 
tify their importance), but we expect that at Tycho’s age, 
X-ray emission is dominated by forward shock transport. 
Given the uncertainty in shock morphology, our model 
seeks only to capture the most relevant physics. More 
sophisticated work may treat, e.g., sphericity, shock pre- 
cursors, anisotr opic diffusion, and in j ection/acceleration 


efficiency (e.g. , Rev ille & Bell 2013 Bykov et al. 2014 
Ferrand et al.|[2U14 and references therein). 


To determine rim profiles and widths, we compute the 
electron distribution using G reen’ s function solut i ons by 
Lerche & Schlickeiser (1980) and Rettig & Pohl (2012), 
with the caveat that D(x)B 2 ( x) is assumed constant; 
we discuss this assumption further below. The solutions 
are fully described in |R14| using notation similar to ours. 
The electron distribution may be integrated over the one- 
particle synchrotron emissivity G{y) to obtain the “to- 
tal” emissivity: 


jv{x) oc / G(y)f(E, x)dE 


( 2 ) 


where y = v / (c\E 2 B) is a scaled synchrotron frequency 
and G{y ) = y f y K 5 / 3 (z)dz with K 5 / 3 (z) a modified 
Bessel function of the second kind (Pacholczyk 1970); 
the constant c\ = 6.27 x 10 18 in CGS units. Integrat- 
ing emissivity over lines of sight for a spherical remnant 
yields intensity as a function of radial coordinate r: 


I u (r) = 2 J j v (r s - \j s 2 T r 2 ) 


ds 


( 3 ) 


where s is the line-of-sight coordinate and r s is shock ra- 
dius. We take the full width at half maximum (FWHM) 
of the resulting intensity profile as our metric for modeled 
rim widths. Using FWHM as opposed to, e.g., full width 
at three-quarters maximum, excludes measured rims and 
model parameters where X-ray intensity does not drop to 
half maximum immediately behind the rim. Our results 
thus focus on the most well-defined rims in Tycho rather 
than the global shock structure. 


2.2. Magnetic fields and damping 

We consider two scenarios for post-shock magnetic 
field: (1) a constant field B(x) = B 0 corresponding to 
loss-limited rims, and ( 2 ) an exponentially damped field 
of form: 

B(x) = ( B 0 - B min ) exp (—x/ a b ) + B min (4) 

following (Pohl et al. |2005 ). Here Bq is the magnetic field 
immediately downstream of the shock, i.e. Bq = B{x = 
0 ) , and ab is an e- folding damping lengthscale; our use of 
B 0 for downstream magnetic field departs from typical 
notation. A typic al lengthscale for ab is 10 16 to 10 17 cm 


(Pohl et al.| 2005 ), corresponding to ~0.1- 1% of Tycho’s 
radius. Hereafter, we report ab in units of shock radius 
r s unless otherwise stated. 

The distinction between damped and loss-limited rims 
is somewhat arbitrary; as Qj, -> oo, model results con- 
verge to loss-limited rims. Moreover, rims much thin- 
ner than the damping length are effectively loss-limited 
as electrons radiate in a nearly constant magnetic field. 
Furthermore, the dependence on observing frequency of 
electron losses and diffusion means that a model may be 
damping limited at one frequency and loss-limited at an- 
other. In the following analysis, we deem all fits with 
finite ab to be damped, but we compare rim widths and 
damping lengths from such fits further below to better 
distinguish damped and loss-limited rim behavior, set by 
a combination of ab, Bq, and other model parameters. 

2.3. Diffusion coefficient 

Most previous work has assumed Bohm-like diffusion 
in plasma downstream of SNR shocks. Bohm diffusion 
assumes that the particle mean free path A is equal to 
the gyroradius r g = E/(eB), yielding diffusion coefficient 
Db = Ac/3 = cE/(3eB ); here c is the speed of light, E 
is particle energy, e is the elementary charge, and B is 
magnetic field. Bohm-like diffusion encapsulates diffu- 
sion scalings of D oc E, introducing a free prefactor 77 
such that A = yr g allows for varying diffusion strength. 
However, Bohm diffusion at 77 = 1 is commonly con- 
sidered a lower limit on the diffusion coefficient at all 
energies. 

We consider a generalized diffusion coefficient with ar- 
bitr ary power law dependence upon energy following, 
e.g., Parizot et al. (2006): 


D(E) = ^ = Vh D B (E h ) 


E 

Eh 


( 5 ) 


where y parameterizes diffusion-energy scaling and 77 now 
has units of erg 1-Al . The right-hand side of equation (l5| 
introduces 77 h, a dimensionless diffusion coefficient scaled 
to the Bohm value at a fiducial particle energy Eh- Note 
that r/h and 77 are related as 77 = rjh(Eh ) 1 ~ fl , and rj = rjh 
for Bohm-like diffusion (/x = 1). For subsequent analy- 
sis, we take fiducial electron energy E 2 = E h correspond- 
ing to a 2 keV synchrotron photon and report results in 
terms of 772 = rjh- Although E 2 varies with magnetic field 
as E 2 oc B ~ 1 / 2 and thus 772 may vary around Tycho’s 
shock for y 7 ^ 1 , tying 77 h to a fixed observation energy 
gives a convenient sense of diffusion strength regardless 
of the underlying electron energies. 

The solutions t o equation (fTl) given by Lerche & 
Schlickeiser (1980) assume D(x)B 2 (x) constant to ren- 
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der equation ([T]) semi-analytically tractable. Although 
this assumption has no obvious physical basis, it con- 
tains the qualitatively correct behavior D constant if B 
is constant and D smaller for larger B. We enforce it by 
modifying th e diffusion coefficient in the damping model 
as, following Rettig & Pohl (2012): 


D(E,x) = 


rjcE 11 

3eBn 


B m i 


Bq — B 
~ Bn 


min g—x/db 


-2 


( 6 ) 


This strengthens the spatial-dependence of the diffusion 
coefficient as compared to the expected D(x ) oc 1/B(x). 

2.4. Electron energy cut-off 

We assume that the DSA process is limited by syn- 
chrotron losses at high energies and hence determine 
the E cut by equating synchrotron loss and diffusive ac- 
celeration timescales. Here we ignore the possibility of 
age-l imited or escape-limited acceleration (e.g.,|Reynolds 
1998). In the former case, low magnetic-field strengths 
could mean that Tycho’s age is less than a synchrotron 
loss time; the maximum energy is obtained by equating 
the remnant age and acceleration timescale. In the lat- 
ter case, the diffusion coefficient upstream may increase 
substantially above some electron energy due to an ab- 
sence of appropriate MHD waves. However, the mag- 
netic field strengths we find below justify the assumption 
of loss-limited acceleration. For low energies and small 
synchrotron losses (cooling time longer than acceleration 
time), electrons are efficiently accelerated; near or above 
the cut-off energy, electrons will radiate or escape too 
rapidly to be accelerated to higher energies and the en- 
ergy spectrum drops off steeply. The cut-off energy is 
given as: 


E cut = (8.3 TeV) 2/(1+M) ( 


Bo 


( 


100 pG 
\ 2/(1+a0 


-i/(i+/d 


10® cm s -1 / 


V 




( 7 ) 


This result is derived by Parizot et al. (2006) for p = 1 
assuming a strong shock with compression ratio 4 and 
isotropic magnetic turbulence both upstream and down- 
stream of the shock. 

The cut-off energy and accelerated electron spectrum 
depend on the ma gnetic f ield, which varies in a damped 


model. Marcowith & Casse (2010) show that various 


turbulent damping mechanisms can modify the acceler- 
ated spectrum, as particles traveling downstream may 
not be effectively reflected back across the shock and 
further accelerated. Nevertheless, we make the simpli- 
fying assumption that particle acceleration is controlled 
by diffusion at the shock and neglect spatially varying 
diffusion and magnetic fields in the acceleration process; 
equation 0 stands as evaluated with shock magnetic 
field strength B 0 . Cut-off energies of ~1 10 TeV can 
be plausibly achieved in the presenc e of d amping due to 
Alfven and magneto-sonic cascades (Marcowith & Casse 
2010 ). 


"As the DSA imposed electron cut-off results in a cut- 
off of SNR synchrotron flux, the synchrotron cut-off fre- 
que ncy v CMt = Cr„Ef,+ B with c m = 1.82x 10 18 in cgs units 
(e.g. Pacholczyk 1970), which is the peak frequency emit- 
ted by electrons of energy E cut , provides an independent 


observable to estimate shock diffusion and is given by: 
z'cut = c m (13.3 erg) 1 ^ (100 pG) (2657 erg 2 ) 1+ " 


10 s cm/s 


0/2 )' 


1 + M 


( 8 ) 


The cut-off frequency is independent of magnetic field 
B for all values of p (but recall that the electron en- 
erg ies associat e d wit h r] 2 will depend on B for / 1). 
Parizot et al. (2006| previously used measurements of 
synchrotron cut-off's to estimate diffusion coefficients in 
Tycho and other historical supernova remnants. 

We point out that the electron spectrum softens down- 
stream of the shock due to synchrotron losses, but the 
local spectrum at any radial position will be a ste eply 
cut-off power law ( Webb et al.|1984 Reynolds|1998 ). No 
steepening from E~ s to FCC+b is observed because the 
exponential cut-off limits the electron spectrum at high 
energies. A homogeneous source of age t in which elec- 
trons are continuously accelerated throughout, with an 
initial straight power law distribution to infinite energy, 
will produce a steepened power law distribution above 
the energy at which the sy nchrotron loss time equals the 
acceleration time (Kardashevj 1962). If we model emis- 
sion without an initial exponential cut-off (i.e. , inject 
assuming a straight power law) in a constant magnetic 
field, then integrated spectra in our model would steepen 
by about one power. But, those assumptions do not ap- 
ply to the current situation of continuous advection of 
electrons, with an energy distribution which is already a 
cut-off power law with expon e ntial cutoff , through 
a reg ion of non-constant magnetic field. See |Reynolds] 
(2009) for a fuller discussion of synchrotron losses in non- 


homogeneous sources. 

2.5. Rim width-energy dependence 

X-ray rim widths are controlled by synchrotron losses, 
particle transport, and magnetic fields immediately 
downstream of the shock, each of which influences rim 
width-energy scaling differently. If diffusion is negligi- 
ble and the downstream magnetic field is constant, loss- 
limited rims narrow with increasing energy as more ener- 
getic electrons radiate and cool more quickly. But diffu- 
sion will dilute this effect: more energetic electrons may 
diffuse further upstream or downstream than would be 
expected from pure advection, smearing out rims and 
weakening energy dependence at higher energies. Mag- 
netic fields damped on a length scale comparable to fila- 
ment widths should also weaken rim width-energy depen- 
dence - if the magnetic field turns off, synchrotron radia- 
tion turns off regardless of electron energy. Additionally, 
once Bq varies, one observes electrons of different energy 
at different distances behind the shock. 

Figure [2] plots model X-ray and radio profiles for a 
range of B 0 and a & values to illustrate how damping and 
magnetic field strength impact width-energy dependence, 
which we will now explore. We discuss and incorporate 
model radio profiles into our analysis in Section [g] 

2.5.1. Undamped models 

Following |R14 we parameterize rim width-energy de- 
pendence in terms of a scaling exponent ttie defined as: 


w(v) oc v 71 


(9) 
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Figure 2. Radio and X-ray rim profiles for a range of magnetic fields ( Bq ) and damping lengths (a. 5 , shaded regions) with fi = 1, 772 = 1, 
and B m in = 5 fiG fixed. X-ray rim energies (0.7, 2, 4 keV) are representative of energy bands used in our rim width measurements. This 
plot summarizes several key features of our model: (1) radio profiles (dashed red) are strongly affected by damping for all values of Bo 
because synchrotron losses are negligible at low electron energies. (2) X-ray profiles (solid green, black, blue) are influenced by synchrotron 
losses even in the presence of strong damping, which can be seen going from left to right (increasing Bo). (3) When X-ray rim widths are 
smaller than a 5 , field damping only weakly affects rim widths and width-energy dependence (top right panels). (4) Strongly damped X-ray 
rims can show significant width-energy dependence in our model due to synchrotron losses beyond a 5 . This occurs for smaller Bo where 
the contrast between Bq and B m i n is less extreme (bottom left panels). 


where w{y) is filament FWHM as a function of observed 
photon frequency u, and the exponent ttte = is 

energy dependent. We may refer to observed photons 
interchangeably by energy or frequency u, but E is re- 
served for electron energy. 

To better intuit the effects of advection and diffu- 
sion on rim widths, we introduce advective and diffusive 
lengthscales for bulk electron transport. These depend 
on electron energy; we write them in terms of the peak 
frequency radiated by electrons of energy E, v = c m E 2 B. 

Ld = v d T synch oc v d B// 3 , 2 v ~ 1/2 (10) 

Zdiff = \J Dr synch oc ?? 1 /2 s -(^+ 5 )/ 4 j / (a*-i)/ 4 (n) 

The characteristic time is the synchrotron cooling time 
Tsynch = 1 /(bB 2 E) with b = 1.57 x 10~ 3 . For fi = 1, Zdiff 

is independent of v and both (diff and Z a d scale as B 0 3 ^ 2 . 
If both diffusion and magnetic Held damping are negli- 
gible and electrons are only loss-limited as they advect 


downstream, toe attains a minimum value tbe — —1/2 
as rim widths are set by l ad > Zdiff- At higher energies 
where Zdiff > Z a d> diffusion increases Toe from —1/2 to a 
value between —1/4 and 1/4 for /.i = 0 and 2 respectively 
(R14 Figure 3). The presence of an electron energy cut- 
off decreases vie slightly in all cases due to the decreased 
number of electrons and hence thinner rims at higher en- 
ergies (R14 Figure 5), but the qualitative behavior is the 
same. 


2.5.2. Field damping effects 

We expect magnetic damping to produce compara- 
tively energy-independent rim widths. If synchrotron 
rim widths are set by magnetic damping at some ob- 
servation energy, then rims will be damped at all lower 
observation energies as well. Then rim widths will be 
relatively constant (small |?71e|) below a threshold en- 
ergy and may decrease, or even increase once advection 
and/or diffusion control rim widths at higher photon en- 


6 


Tran et al. (September 2, 2015) 


ergies (advection: Z ac j < o&; diffusion: Zdiff > Gd,Ob). 
Thus, we intuit that rim widths should roughly scale 
as w ~ min (af,, max (/ a d, ?diff))- This is correct except 
for one key region of parameter space: strong damp- 
ing with weak magnetic field amplification, where syn- 
chrotron losses downstream of the FWHM create energy 
dependent widths even when cif, <C Cardiff- 
This counter-intuitive energy dependence for strongly 
damped models occurs when rim brightness remains 
above half-maximum within ~a& of the shock. Farther 
downstream, synchrotron losses in the reduced magnetic 
field ~B m drive intensity down to half-maximum; losses 
in a nearly constant field cause energy dependence de- 
spite magnetic damping at the shock. Synchrotron losses 
far from the shock will only affect the spectrum at higher 
frequencies; at lower frequencies, the spectrum will be 
effectively constant after the magnetic field’s initial de- 
cay. Thus, there will be a characteristic frequency be- 
low which the FWF1M ceases to be defined. The energy 
dependence of rim FWFIMs becomes large and diverges 
near this frequency. The characteristic frequency and 
value of we are sensitive to our definition of rim width 
(FWHM), but the physical behavior we describe (diver- 
gence of Toe at characteristic frequency, for appropriate 
damping parameters) will occur regardless of our defini- 
tion of rim width. 


3. OBSERVATIONS 
3.1. Data and region selections 

We measured synchrotron rim full widths at half max- 
imum (FWHMs) from an archival Chandra ACIS-I ob- 
servation of Tycho (RA: 00 h 25 m 19!0, dec: +64°08 , 10'.'0; 
J2000) between 2009 Apr 11 and 2009 May 5 (PI: J. 
Hugh es: ObsIDs: 10093-10097, 10902-10906); |Eriksen 
et al. (2011) present additional observation information. 
The total exposure time was 734 ks. Level 1 Chan- 
dra data were reprocessed with CIAO 4.6 and CALDB 
4. 6. 1.1 and kept unbinned with ACIS spatial resolution 
0.492". Merged and corrected events were divided into 
five energy bands: 0.7-1 keV, 1- 1.7 keV, 2-3 keV, 3-4.5 
keV, and 4.5-7 keV. We excluded the 1.7-2 keV energy 
range to avoid Si XIII (Hea) emission prevalent in the 
remnant’s thermal ejecta which might contaminate our 
nonthermal profile measurements. 

We selected 20 regions for profile extraction around 
Tycho’s shock (Figure [I]) based on the following criteria: 
(1) filaments should be clear of spatial plumes of thermal 
ejecta in Chandra images, which rules out, e.g., areas of 
strong thermal emission on Tycho’s eastern limb; (2) fil- 
aments should be singular and localized, so multiple fila- 
ments should either not overlap or completely overlap; 
(3) filament peaks should be evident above the back- 
ground signal or downstream thermal emission (rules 
out faint southern filaments). We accepted several re- 
gions with poor quality peaks in the lowest energy band 
(0.7 1 keV) so long as peaks in all higher energy bands 
were clear and well- fit. We grouped regions into 5 fila- 
ments by visual inspection of the remnant. Within each 
filament, we chose region widths to obtain comparable 
counts at the thin rim peak. All measured rim widths are 
at least 1". The narrowest rim widths may be slightly 
over-estimated due to the Chandra point-spread func- 
tion (PSF) at ~4' off the optical axis, which has FWHM 


< 1.4" at 6.4 keV and < 1" at 1.5 keV. For simplicity, 
we neglect PSF effects in our analysis. Our observations 
of rim width-energy dependence are thus somewhat con- 
servative at the highest energy 


3.2. Filament spectra 


We extracted spectra at and immediately behind thin 
rims in each region ( “rim” , “downstream” spectra respec- 
tively) to confirm that rim width measurements are not 
contaminated by thermal line emission. The two extrac- 
tion regions are deter min ed by our empirical fits of rim 
profile shape (Section |3.3[ ). The rim section is the small- 
est sub-region containing the measured FWHM bounds 
from all energy bands. The downstream section extends 
from the interior thin rim FWHM to the intensity mini- 
mum behind the rim (specifically, the do wnst ream profile 
fit domain bound described in Section 3.3). To illus- 
trate our selections, Figure [3] plots example rim profiles 
(4.5-7 keV) with the downstream and rim sections high- 
lighted. 

Spectra were binned to a minimum of 15 cts/bin. We 
extracted background spectra from circular regions (ra- 
dius ~30") around the remnant’s exterior; for each re- 
gion’s rim and downstream spectra, we subtracted the 
closest background region’s spectrum. 

We fit each region’s rim and downstream spectra 
to an absorbed power law model (XSPEC 12.8.1, 
phabs*powerlaw) between 0.5-7 keV with photon index 
r, hydrogen column density _/Vh, and a normalization as 
free parameters. Table [l] lists best fit parameters and re- 
duced x 2 values for all regions. Rim spectra are well-fit 
by the power law model alone; the best fit photon in- 
dices (2.4-3) and column densities (0.6-0. 8 x 10 22 cm -2 ) 
are consistent with previous spect ral fits to Tych o’s non- 


thermal rims (Hwang et al. 2002| Cassam-Chenai' et al. 


2007). 

Downstream spectra are poorly fit by the absorbed 
power law model due to thermal contamination from 
Si XIII and S XV Hea line emission at 1.85 and 2.45 
keV. To confirm that thermal emission is dominated by 
these two lines near the shock, we also performed fits with 
(1) both lines excised (1. 7-2.0 keV, 2. 3-2. 6 keV counts 
removed) and (2) with both lines fitted to Gaussian pro- 


files. Fits with lines excised yield x? e d values between 
1-5. Fits with lines fitted to Gaussian profiles yield Xred 
values 0.83-1.6. In both fits (lines excised or modeled), 
we find somewhat smaller best fit column densities (0.3- 
0.8 x 10 22 cm” 1 ) but similar best fit photon indices (2.6- 
3.1), compared to those of the rim spectra. The consis- 
tent photon indices indicate that the same synchrotron 
continuum is present beneath thermal li ne em ission. 

We also fitted “rim” spectra (Section |3.2[ ) to the ab- 
sorbed XSPEC model srcut, modified to fit in log- 
frequency space, srcut models a power law X-ray syn- 
chrotron spectrum set by a radio spectral index a with 



lists best spectrum fit parameters for each region. The 
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Figure 3. Spectra and fits from Regions 1 (top) and 16 (bottom) show varying rim morphology; Region 1 shows a rim where the 0.7-1 keV 
peak could not be fit. Left: 4.5-7 keV profiles with downstream (blue) and rim (grey) sections highlighted. Intensity is in arbitrary units 
(a.u.). Middle: downstream spectra with absorbed power law fit; Si and S lines at 1.85, 2.45 keV are clearly visible. Right: rim spectra 
with absorbed power law fit show that rims in each region are likely free of thermal line emission. 


Table 1 

Absorbed power law spectrum fit parameters 


Downstream, power-law Rim, power-law Rim, srcut 


Region 

JV H 

r 

Xred ( d ° f ) 

JV H 

r 

X?ed (dof) 

JVh 

V cut 

Xred ( d ° f ) 


(10 22 cm — 2 ) 

(-) 


(10 22 cm — 2 ) 

(-) 


(10 22 cm -2 ) 

(keV /h) 


1 

0.53 

2.72 

2.25 (186) 

0.72 

2.84 

1.20 (284) 

0.62 

0.30 

1.17 (284) 

2 

0.68 

2.99 

5.07 (178) 

0.69 

2.83 

1.10 (202) 

0.60 

0.30 

1.08 (202) 

3 

0.67 

2.96 

1.97 (186) 

0.77 

2.80 

1.15 (167) 

0.68 

0.33 

1.14 (167) 

4 

0.59 

2.97 

1.43 (163) 

0.70 

2.84 

1.21 (278) 

0.61 

0.29 

1.15 (278) 

5 

0.62 

2.93 

4.76 (265) 

0.73 

2.88 

1.18 (255) 

0.64 

0.27 

1.11 (255) 

6 

0.68 

3.00 

2.12 (200) 

0.74 

2.85 

0.96 (231) 

0.64 

0.29 

0.93 (231) 

7 

0.65 

3.02 

1.00 (142) 

0.82 

2.97 

1.14 (224) 

0.71 

0.23 

1.14 (224) 

8 

0.74 

2.93 

1.38 (170) 

0.75 

2.71 

0.98 (198) 

0.66 

0.41 

0.96 (198) 

9 

0.78 

3.03 

1.12 (157) 

0.82 

2.83 

0.90 (175) 

0.73 

0.30 

0.88 (175) 

10 

0.62 

2.86 

1.40 (220) 

0.77 

2.76 

0.98 (164) 

0.68 

0.36 

0.96 (164) 

11 

0.67 

2.94 

2.56 (137) 

0.69 

2.60 

1.10 (153) 

0.61 

0.55 

1.07 (153) 

12 

0.61 

2.79 

2.65 (137) 

0.64 

2.44 

0.90 (172) 

0.57 

0.88 

0.90 (172) 

13 

0.61 

2.98 

3.12 (198) 

0.67 

2.73 

1.12 (235) 

0.59 

0.38 

1.09 (235) 

14 

0.46 

2.93 

1.37 (148) 

0.63 

2.93 

0.96 (167) 

0.54 

0.23 

0.95 (167) 

15 

0.43 

2.92 

1.33 (150) 

0.65 

2.84 

1.05 (183) 

0.57 

0.29 

1.03 (183) 

16 

0.48 

2.94 

2.04 (189) 

0.67 

2.80 

1.13 (182) 

0.58 

0.32 

1.12 (182) 

17 

0.48 

2.86 

1.70 (188) 

0.68 

2.83 

0.97 (187) 

0.59 

0.30 

0.94 (187) 

18 

0.44 

2.87 

1.91 (200) 

0.64 

3.02 

1.20 (220) 

0.55 

0.19 

1.13 (220) 

19 

0.40 

2.84 

1.31 (133) 

0.66 

2.78 

1.01 (157) 

0.57 

0.34 

0.99 (157) 

20 

0.40 

2.75 

3.01 (140) 

0.63 

2.81 

1.11 (192) 

0.55 

0.31 

1.07 (192) 


Note. — Absorbed power law fit parameters are photon index T and hydrogen column density Nu. srcut fits performed in log-frequency space; 
h is Planck’s constant and v cut a cut-off frequency. Horizontal rules group individual regions into filaments. 
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fitted cut-off fre quency is typic ally 0.3 keV /h, consistent 
with fits by Hwang et al. (2002), but Regions 11, 12 have 
unusually high cut-off frequencies 0.55, 0.88 keV/h con- 
sistent with harder rim spectra. 

Our spectral fitting confirms that all selected region 
are practically free of thermal line emission, as already 
suggested by visual inspection (Figure[l]). Excluding 1.7 
2 keV photons in rim width measurements further lim- 
its thermal contamination as 1.85 keV Si line emission 
is over a third of Tycho’s th ermal flux as detected by 
Chandra (Hwang et al. 2002 1 . 


3.3. Filament width measurements 

We obtained radial intensity profiles in five energy 
bands from ~10-20" behind the shock to ~5-10" in front 
for each region. To increase signal-to-noise, we integrate 
along the shock (5-23") in each region. Plotted and 
fitted profiles are reported in vignetting and exposure- 
corrected intensity units; error bars were computed from 
raw counts assuming Poisson statistics. Intensity profiles 
peak sharply within ~2-3" behind the shock, demarcat- 
ing the thin rims, then fall off gradually until thermal 
emission picks up further behind the shock. 

We fitted rim profiles to a piecewise two-exponential 
model: 

h(r) = j A “ eXP ( “«* ) +Cu ’ (12) 

[Ajexp^j +C d , r < r 0 

where h(r ) is profile height and r is radial distance from 
remnant center. The rim model has 6 free parameters: 
A-u , Cq , Uj u ■ w'd , Cu i and C^, A u ( C u Cfi) en- 

fo rces continuity at r = rp. Our model is similar to that 
of lBam ba et al. ( 2003 2005 1 and differs slightly from that 
of |R14 To fit only the nonthermal rim in each intensity 
profile, we selected the fit domain for each profile as fol- 
lows. The downstream bound was set at the first local 
data minimum downstream of the rim peak, identified by 
smoothing the profiles with a 21-point (~10") Hanning 
window. The upstream bound was set at the profile’s 
outer edge. Figure [4] illustrates the fit domain selections 
for two example regions. 

From the fitted profiles we extracted a full width at 
half maximum (FWHM) for each region and each en- 
ergy band after subtracting a constant background term 
min (C u ,Cd)- We could not resolve a FWHM in 8 of 20 
regions at 0.7 1 keV (Table [2); in these regions, either 
the downstream FWHM bound would extend outside the 
fit d oma in or we could not find an acceptable fit to equa- 
tion ( |T2| ) . We were able to resolve FWHMs for all regions 
at higlfer energy bands (1-7 keV). 

To estimate FWHM uncertainties, we horizontally 
stretched each best fit profile by mapping radial coor- 
dinate r to r'[r) = r( 1 + £(r — r o )/(50" — r 0 )) with 
£ an arbitrary stretch par ame ter and ro the best fit 
rim center from equation (12), yielding a new profile 
h' (r) = h(r'(r)). We varied £""{and hence rim FWHM) to 
vary each profile fit y 2 by 2.7 and took stretched FWHMs 
as upper/lower bounds on reported FWHMs. This pro- 
cedure again follows |Rl4| 

3.4. Filament model fitting 


We fit model FWHM predictions given by equation Q 
to measured rim widths as a function of energy by vary- 
ing several physical parameters: magnetic field strength 
B 0 , normalized diffusion coefficient 772 , diffusion-energy 
scaling exponent /z, and minimum field strength B m ; n 
and lengthscale a& for a damped magnetic field. We 
mapped each width measurement to the lower energy 
limit of its energy band; e.g., 0.7 1 keV is assigned to 
0.7 keV and fitted to model profile widths at 0.7 keV. 
Width errors in our least squares fits average the positive 
and negative errors on each FWHM measurement. For a 
given set of model parameters, we numerically computed 
intensit y pr ofiles and hence model FWHMs as detailed in 
Section [2d] we then used a Levenberg-Marquardt fitter 
to seek model parameters yielding best fit FWHMs. To 
assist the nonlinear fitting, we tabulated model FWHM 
values on a large grid of model parameters and used best 
fit grid parameters as initial guesses for fitting. We re- 
quired 772 to be positive and deemed best fit values with 
772 > 10 5 and Bq > 10 rnG to be effectively uncon- 
strained. In subsequent analysis we focus on fits with 
/z = 1 and 772 = 1 fixed, though we discuss the effects of 
varying both /z and 772 (and fitting for 772 ) as well (recall 
that for /z = 1 , 77 is energy- independent). 

The purely loss-limited model (constant downstream 
magnetic field Bq) has three parameters /z, 772 , and Bq. 
To make nonlinear fitting tractable, we fixed /z in all 
fits and considered /z = 0, 1/3, 1/2, 1, 1.5, and 2. In 
particular, nonlinear diffusion-energy scalings with 77 = 
1/3 and 1/2 may arise from Kolmogor ov and K raichn an 
turbulent energy spectra respectively ( Reynolds| [2004 ) . 

For damped magnetic field fits, we held the remnant 
interior field strength B min constant at 5 ^zG, sl ightly 
higher than typical intergalactic values of ~2-3 /zG (|Lyne 
& Smith|| 1989 |Han et a.l.||2006[). We stepped at, through 
14 different values between (J.5 and 0.002 (sampling most 
finely between 0.01 and 0.002). To ensure that damped 
fits generate rims influenced by magnetic damping, we ar- 
bitrarily require that the rim FWHM at 2 keV be strictly 
greater than the fit value of Of,; we revisit this require- 
ment below. For best fits, we report the value of ab 
yielding the smallest y 2 value, with the caveat that that 
our ab sampling is relatively coarse. 

Predicted rim widths are subject to resolution error 
in the numerical integrals (discretization over radial co- 
ordinate, line-of-sight coordinate, electron distribution, 
Green’s function integrals). We chose integration reso- 
lutions such that the fractional error in model FWHMs 
associated with halving or doubling each integration res- 
olution is less than 1% for the parameter space relevant to 
our filaments. The maximum resolution errors in a sam- 
ple of parameter space are typically 0.1-1%, but mean 
and median errors are typically an order of magnitude 
smaller than maximum errors. 


4. RESULTS 
4.1. Rim widths 

Measured rim widths decrease with energy in most re- 
gions and energy bands. Table [2] reports FWHM mea- 
surements for all of our regions. We also report ttie values 
for all but the lowest energy band, computed point-to- 
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Figure 4. Best fit profiles with measured FWHMs demarcated for each energy band in Region 1 (top) and Region 16 (bottom). Hereafter, 
we use Regions 1 and 16 to illustrate results for profiles of differing quality and absolute width. We could not measure a 0.7-1 keV FWHM 
in Region 1, reflected in Table [2] Data in red were excluded from profile fit domains as described in text. 


point between discrete energy bands as: 


m E (E 2 ) 


hx(w 2 /w\) 
In {E 2 /E 1 ) 


(13) 


where w\,w 2 and Ei,E 2 are FWHMs and lower energy 
values for each energy band - e.g., m E at 1 keV is com- 
puted using FWHMs from 0.7 1 keV and 1 1.7 keV, 
with E 2 = 1 keV and Ei = 0.7 keV. Errors on m E are 
propagated in quadrature from adjacent FWHM mea- 
surements. 

Although the measurement scatter is quite large, 
shown dramatically in the point-wise computed m E val- 
ues, the mean rim width decreases consistently with in- 
creasing energy. Furthermore, mean m E values are con- 
sistently negative and tend smoothly towards 0 (weaker 
energy-dependence) with increasing energy. Errors on 
FWHM measurements are typically < 10%, reflecting 
the high quality of the underlying Chandra data. Scat- 
ter in FWHM measurements may be attributed in part 
to ( 1 ) our measurement procedure, which depends on an 
empirical choice of profile fit function, and ( 2 ) variation 
in Tycho’s rim morphology (e.g., Figure|4|. 

4.2. Model fit results 

Best fit results for loss-limited and damped cases are 
summarized in Figure [5] and Table [3] with fx = 1 and 
772 = 1 both fixed. In all tables, we report values of y 2 
rather than y 2 ed to ease comparison of different fits, as 
we often manually stepped parameters that were held 
constant in fitting. 

Table [4] reports model fits with r ] 2 fixed to values de- 
rived from srcut spectrum fits for z/ cut and hence r ] 2 
from equation (18]), with Bq varying freely; cutoff-derived 
772 values of ~lC) require increased magnetic fields com- 
pared to the 772 = 1 case if rims are purely loss-limited. 
We derive different 772 values for varying 77 and find that 
fitted B 0 varies only weakly with 71 ; for example, Re- 
gions 1 and 16 have best fit values of Bq increasing over 
252-266 77 G and 660-700 77 G respectively as 77 ranges 
between 0 - 2 , for loss-limited fits. 

Magnetically damped rims are able to fit width-energy 
dependence in our measurements at least as well as 
purely loss-limited rims for an acceptable range of ab 
values. As loss-limited models are a subset of damping 


models, this is expected. But fits with damping lengths 
small enough to exert influence on rim widths are also 
permissible. In both loss-limited and damped fits, we do 
not observe systematic variation of fit parameters with 
azimuth around the remnant. 

4.3. Fitting parameter degeneracy 

We briefly explore the effects of varying 772 and other fit 
parameters, before focusing on results with 77 = 772 = 1 
for simplicity. Fixing 772 from measured ix cu t values does 
not yield any obvious insight. But, diffusion coefficients 
computed in this manner may give more credible derived 
estimates of Bq if the DSA assumptions invoked are ac- 
curate. 

Loss-limited model values of Bq and rj 2 may covary 
without strongly altering fit quality in many regions. If 
diffusion is the primary control on rim width i.e. , > 

lad ~ then 772 and B 0 become degenerate as the product 


-3/2 


rj 1 / 2 B i 

(equation (111). If advection is the primary control on 


exerts most control on rim width w ~ Zdiff 


rim width (widths narrow rapidly with energy; i.e., m E 
—0.5), then 772 1 becomes unimportant and fits are 

well-behaved with effectively one free parameter. 

We also performed loss-limited model fits with 77 fixed 
between 0 and 2 and both rj 2 and Bq free. Fits with 
77 < 1 generally yield larger parameter values and errors 
for both 772 and Bq, and they are more likely to be ill- 
constrained (e.g., 772 > 10 3 or B 0 > 10 3 77G). Fits to the 
same data with varying 77 can yield rj 2 varying by 1 -2 or- 
ders of magnitude. The % 2 values for individual fits are 
variable and large (3> 1 ), so we cannot favor or disfavor 
particular values of 77 and 772. But, neglecting the mag- 
nitude of our x 2 values, values of 77 > 1 are qualitatively 
favored by y 2 in most regions. This trend may be par- 
tially an artifact of the correlation between Bq and rj 2 as 
they are not entirely independent parameters, but they 
are less correlated for larger 77; fits at smaller 77 may have 
fewer (non-integer) degrees of freedom, partially offset- 
ting our observations. 

Damping modifies the degeneracy in B 0 and r] 2 . For 
small values of ab (strong damping), increased diffu- 
sion ( 772 ) can cause rim widths at all energies to narrow 
counterintuitively. Speculatively, spatial variation of the 
diffusion coefficient may oppose downstream advection 
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Table 2 

Measured full widths at half max (FWHMs) for all regions. 





FWHM (arcsec 

) 




m E ( 

» 


Region 

Band 1 

Band 2 

Band 3 

Band 4 

Band 5 

Bands 1-2 

Bands 2—3 

Bands 3—4 

Bands 4—5 


(0.7-1 keV) 

(1-1.7 keV) 

(2-3 keV) 

(3-4.5 keV) 

(4.5-7 keV) 

(1 keV) 

(2 keV) 

(3 keV) 

(4.5 keV) 

1 


8-801“;;® 

6 34 + 0 - 26 
U - Oit _ 0.21 

7.40l“ o ; 3 2 “ 

5 57 ”*"°' 47 

-0.42 


-0.47 ± 

0.06 

0.38 ± 0.13 

-0.70 ± 0.22 

2 


4-22i ° 0 ,; 2 9 

2 36~*” 012 

z.oo _ 0 09 

3.001“,;;“ 

4 n+0-34 
J " L — 0.30 


-0.84 ± 0.08 

0.59 ± 0.16 

0.77 ± 0.23 

3 


2 47 + 0 - 08 
Z -^'-0.07 

1 70 + 0.09 
1-7o_ 0 .07 

2 10" 1 " 011 
z - ±u -o.n 

1 q 9 + 0 . 10 

1 - OZ -0.09 


-0.47 ± 

0.08 

0.41 ± 0.17 

-1.15 ± 0.22 

4 

5 85"*" 0 ' 37 

O.OO _0 33 

4 S5 "*" 0 09 
^.oo_o os 

3-26l“ 0 ;“ 

3 59+° 12 
o.uy _ 0 11 

3 20 ”*”°' 21 
— 0.18 

-0.83 ± 0.18 

-0.41 ± 

0.05 

0.31 ± 0.11 

-0.35 ± 0.17 

5 


4 52 -1 ” 0 ' 11 
^• OZ - 0.12 

3-o6is.;; 

q 9 C+O .15 
°- zo -0.13 

3 04 +o ' 21 
°- u ^_0.18 


-0.56 ± 

0.06 

0.15 ± 0.14 

-0.17 ± 0.19 

6 

2 48" 1 " 0 ' 18 
^•^0-0.18 

9 09 + 0.05 
z - oz - 0.06 

2 98" 1 " 011 
^.yo _ 0 09 

2 n c s+ 0 08 
z.uo _ 0 09 

9 91 +0.15 
^•^-L_o.14 

-0.19 ± 0.21 

0.36 ± 

0.06 

-0.92 ± 0.13 

0.18 ± 0.19 

7 

2 69" 1 " 0 ' 20 
z.oy _ 0 i 7 

9 qq + 0.05 

Z.OO _ 0 05 

9 91 +0.08 
z ' ol -0.08 

1 01 +0.09 
l-ol_ 0.07 

1 oq+ 0-11 

l.OO-Q 08 

-0.39 ± 0.20 

-0.01 ± 

0.06 

-0.60 ± 0.14 

0.02 ± 0.17 

8 

0 oo"l" 0'21 

Z.OO — 0 20 

0 »7q+0.08 
1 -0.08 

9 qo+0.10 
z.oo _ 0 09 

2 lO" 1 ’ 0,10 

_LU — 0.09 

9 q7 + 0 . 20 

z ot -0.17 

0.43 ± 0.26 

-0.19 ± 

0.07 

-0.30 ± 0.15 

0.29 ± 0.22 

9 

2 16+ 0 - 24 
z - lo - 0.23 

9 oc+0.07 

Z.OO _ 0 06 

2 47+ 011 
- 0.11 

1 9 I + 009 

± — 0.09 

2 20 -1-0 ' 17 
Z ' ZU -0.16 

0.24 ± 0.31 

0.07 ± 

0.07 

-0.63 ± 0.16 

0.34 ± 0.22 

10 

o oo-|-0.24 
Z.OO_o 23 

1 gg+0.07 

T.yy _ 0 06 

1 70+0.09 
' u -0.08 

1 59+ 009 
T.oy_ 0 os 

1 ^ 8 +013 
i 0 °- 0.12 

-0.50 ± 0.29 

-0.18 ± 

0.08 

-0.24 ± 0.18 

-0.02 ± 0.23 

11 


0 90 + 0-15 
— 0. 13 

9 C 9 +O.I 6 

z - oz -0.13 

1 90 +0 14 

i.yu_o 13 

3 09"*~°‘ 45 
o.uy _ 0 38 


-0.36 ± 

0.10 

-0.70 ± 0.22 

1.21 ± 0.37 

12 


3-861“;" 

2 61+ 015 
z,ol -0.13 

q r» 9 + 0.22 
^• u ^_ 0.21 

9 9q+0.21 


-0.56 ± 

0.10 

0.36 ± 0.22 

-0.74 ± 0.27 

13 

2-85 ±S: a i? 

2-431“““ 

9 qc+0.08 
^• ciD -0.05 

1 95+ 009 
T.yu _ 0 10 

1 84+ 0 - 11 
J --°^-0.14 

-0.45 ± 0.20 

-0.04 ± 

0.05 

-0.47 ± 0.13 

-0.15 ± 0.20 

14 

2-86 t“" 

2 42+ 0 06 

-0.04 

9 qo+0.08 

z -^-0.07 

9 qO+0.10 

Z.OO _ 0 os 

2 19"*"°' 12 
^•-L^-O.IO 

-0.47 ± 0.17 

- 0.12 ± 

0.06 

0.17 ± 0.12 

-0.20 ± 0.15 

15 

9 '71 +0. 17 
1 — 0. 16 

1 99 + 0 05 

J-.£7C7_o 04 

1 . 8 U _ 0 05 

00 

1 + 

0 0 
b b 

01 

1 C 9 + 0.09 
- L - oz _0.08 

-0.85 ± 0.18 

-0.15 ± 

0.05 

0.09 ± 0.11 

-0.51 ± 0.16 

16 

1 87+0.14 
—0.13 

1 70 + 0.04 
'°-0.03 

1 EI9+0.06 

1 - oz -0.05 

1 9^+ 9 - 98 
1 - zo -0.04 

1 9 q "+ 0.08 
1ZO -0.06 

- 0.22 ± 0.21 

-0.18 ± 

0.06 

-0.49 ± 0.13 

-0.04 ± 0.17 

17 

1 6 C S + 0 - 13 
i.oo _ 0 12 

1 92" 1 " 0 ' 05 
0.05 

1 54+ 006 

1 -°^-0.07 

1 45 + 0 07 
1 -^°-0.06 

2 05~*”° 16 

z - uo -0.14 

0.43 ± 0.22 

-0.31 ± 

0.07 

-0.16 ± 0.15 

0.86 ± 0.21 

18 


4 45+0 - 13 

q 1 q+0. 17 
-0.16 

2 96~*~°‘ 20 
z.yo _ 0 19 

1 ck+ 0-21 
l.DO-o 16 


-0.49 ± 

0.09 

-0.17 ± 0.21 

-1.45 ± 0.32 

19 


9 ofi+0.08 

Z.OU — o 06 

0.08 

2 16" 1 " 0 12 
z - lo - 0.11 

1 60" 1 " 017 
l.OU-o 14 


-0.02 ± 

0.08 

-0.13 ± 0.17 

-0.74 ± 0.27 

20 

4 81 +0 - 31 
^• O1 _0.31 

1 84" 1 " 0 ' 06 

-L-o^_0.03 

1 q7 - I - 0.08 
i -°' -0.06 

1 ^5+0-07 

i.oo_o 06 

2 14 +0 ’ 23 
^•-*■^-0.23 

-2.68 ± 0.19 

0.02 ± 

0.07 

-0.44 ± 0.14 

0.77 ± 0.28 

Mean 

2.89 ± 0.35 

3.11 ± 0.37 

2.53 ± 0.23 

2.47 ± 0.30 

2.35 ± 0.23 

-0.46 ± 0.24 

-0.25 ± 

0.06 

-0.14 ± 0.10 

-0.09 ± 0.15 

Note. 

— Mean values computed for all regions; 

mean ra.E values are averages 

for region ttie 

values (i.< 

2 ., not 

computed from 

mean FWHMs) 


Errors on mean values are standard errors of the mean. Horizontal rules group individual regions into filaments. 


Table 3 

Best width-energy fit parameters, 

H S3 772 = 1 

Table 4 

Best model fits for all regions, 772 derived from srcut fits, (i 

Region 

Loss-limited 

Damped 


Regioi 


Loss-limited 

Damped 


B 0 (l*G) 

X 2 B 0 ( M G) 

x 2 

&b 

1 rn B 0 (/xG) 

x 2 Bo (A»G) 

x 2 

CLb 

1 

182 

35.1 

28 

24.3 

0.008 

1 

12.0 

256 

29.1 

19 

25.5 

0.008 

2 

312 

80.9 

25 

74.4 

0.003 

2 

11.8 

443 

99.4 

18 

72.0 

0.003 

3 

426 

21.5 

27 

23.9 

0.002 

3 

10.6 

599 

31.8 

19 

29.3 

0.002 

4 

284 

50.9 

28 

29.9 

0.004 

4 

12.0 

402 

36.2 

19 

34.5 

0.004 

5 

288 

19.3 

25 

14.5 

0.003 

5 

12.9 

419 

32.5 

19 

19.1 

0.004 

6 

410 

173.7 

68 

65.2 

0.004 

6 

11.9 

568 

89.7 

147 

61.9 

0.006 

7 

418 

50.9 

291 

11.6 

0.006 

7 

15.3 

617 

11.7 

473 

11.1 

0.009 

8 

388 

47.7 

138 

10.0 

0.005 

8 

8.4 

511 

15.1 

333 

10.0 

0.008 

9 

414 

65.1 

78 

14.2 

0.004 

9 

11.2 

573 

23.3 

66 

13.7 

0.005 

10 

466 

17.1 

29 

1.5 

0.002 

10 

9.4 

631 

1.4 

478 

2.0 

0.007 

11 

355 

12.0 

317 

15.6 

0.010 

11 

6.1 

452 

18.1 

21 

17.5 

0.003 

12 

317 

9.6 

26 

10.4 

0.003 

12 

3.7 

375 

12.6 

21 

9.7 

0.003 

13 

400 

52.7 

258 

10.4 

0.006 

13 

8.4 

529 

13.8 

368 

9.9 

0.008 

14 

383 

102.9 

60 

12.7 

0.004 

14 

11.9 

547 

28.6 

120 

12.4 

0.006 

15 

431 

70.3 

48 

20.0 

0.003 

15 

9.5 

598 

23.1 

426 

20.0 

0.007 

16 

493 

15.0 

434 

4.0 

0.006 

16 

8.4 

670 

6.3 

22 

5.8 

0.002 

17 

467 

61.6 

232 

29.0 

0.004 

17 

9.0 

637 

32.6 

32 

29.4 

0.003 

18 

283 

28.1 

24 

27.5 

0.003 

18 

14.6 

432 

72.4 

17 

16.6 

0.003 

19 

401 

36.9 

78 

8.3 

0.004 

19 

8.9 

542 

12.8 

209 

7.7 

0.006 

20 

463 

121.8 

58 

93.8 

0.003 

20 

9.6 

631 

96.0 

29 

94.2 

0.003 


Note. — Fits for Regions 1-3, 5, 11, 12, 18, and 19 have 3 degrees Note. — 772 values are computed from equation 0 and held fixed 

of freedom; all others have 4. The choice of a best a b value may be in model fits. All comments for Table [3] apply to srcut fits as well, 

construed as removing one additional dof. Damped fits require a b to 
be smaller than the FWHM at 2 keV in order to rule out effectively 
loss-limited fits with large a 5 . 
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Figure 5. Rim width fits as a function of energy for loss-limited and damped fits with /i = 1 and 772 = 1 fixed for all regions, from best 
fit parameters in Table |3[ Damped model predictions (e.g., Region 5) are not given at low energies if FWHMs cannot be calculated for 
model profiles (i.e., modeled intensity behind thin rim exceeds half-maximum of rim peak within model domain of ~20 // ). Ordinate ( y ) 
axis limits vary between subplots and are offset from the origin to better show model predictions and data variation. 


through a decreased effective velocity Vd — ( dD/dx ) in 
our transport equation: 


Vd - 


dD(x) 


dx 


= K 0 E- s e- E/E ^S{x), ( 14 ) 


requiring thinner rims as 77 2 and hence dD(x) / dx in- 
crease. Moreover, the effective velocity may impact rim 
widths even if advection dominated. In practice, if 772 


varies freely in damped fits, we observe that smaller val- 
ues of Of, permit and favor smaller best fit values of both 
B 0 and 772. 

Best fit Bq values are smallest for 772 approaching 0 
and /i = 2 ; intuitively, /1 > 1 strengthens diffusion at en- 
ergies > 2 keV for fixed 772, permitting a smaller best fit 
772 and hence smaller B 0 . Fixing 772 = 1 (Table pll), how- 
ever, ties the range of B 0 values to the range of observed 
rim widths. Our minimum loss-limited values of Bq 
are consistent with prior estimates of ~ 200-300 fiG for 
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Figure 6. Model predictions illustrate weak-field (Region 1) and 
strong-field damping (Region 16) for damped best fit parameters 
with fi = 1 and p 2 = 1 fixed. In X-ray energies (0.7— 4.5 keV) model 
profiles are not strongly energy dependent, but weak-field damping 
profiles evolve and show no measurable FWHM at sufficiently low 
energies. Profiles are normalized to peak thin rim intensity, and 
shaded regions indicate damping lengthscale at, . Model parameters 
are given in Table [3] 


advection-dominated transport (I Volk et al. 

2005 

Pari- 

zot et al. 2006 

Morlino & Caprioii 2012). It loss-limited, 


Tycho's rims require strong magnetic held amplification 
to '"TOOx typical galactic field values of ~2-3 /zG, versus 
the expected 4x amplification from a strong shock. 


Magnetic damping fits permit much smaller values of 
B 0 , as expected. The minimum value of B 0 is around 
20 /iG, which would require no field amplification beyond 
that from a strong shock. Fixing B m ; n = 2 p , G instead 
of 5 /zG permits smaller fit B 0 values, though fit values 
for both 772 and B 0 still display considerable scatter. 

5. DISCUSSION 

5.1. Damping is poorly-constrained, by X-ray 
width-energy dependence 

Width-energy dependence is not as sensitive a discrim- 
inant between damping and loss-limited models as may 
be intuitively expected. Both loss-limited and damped 
fits do not perfectly capture sharp drop-offs (especially 
between 0.7 1 keV and 1-2 keV), and the data show 
large scatter (Figure |5j. In most regions, best fit width- 
energy curves have toe ~ —0.2. Only sharp decreases 
in rim width (toe ~ —0.5) can disfavor damping as a 
control on rim widths (e.g., Region 18). 

When width-energy dependence is weak, both loss- 
limited and damped models give best fits with similar 


profiles - amplified magnetic fields in both models cause 
rim intensities to drop sharply. When width-energy de- 
pendence is stronger, damping profiles yield energy de- 
pendent rims if magnetic fields are small (say, < 50 /zG) 
and the damping lengthscale is much smaller than rim 
FWHMs. Figure [6] shows this contrasting behavior using 
best fit parameters for Regions 1 and 16. 

The best damped fit parameters for Regions 1-5, 10, 
12, and 18 predict a 1.375 GHz radio profile that does 
not drop below 50% of peak intensity within ~20" down- 
stream of the shock, which we treat as having no mea- 
surable FWHM. These regions are all best fit with field 
B 0 < 40 /zG. We refer to this model behavior as 
“weak-field” damping, associated with weak magnetic 
fields and stronger emission intensity immediately down- 
stream of the thin rim. As observation frequency de- 
creases, rim and trough contrast decreases, causing rim 
FWHMs to increase and eventually become unmeasur- 
able. Rim width-energy dependence strengthens dra- 
matically (toe — t — oo), permitting model fits to repli- 
cate strong width-energy dependence in observed rim 
widths. Weak-field fit parameters require energy losses 
downstream of the shock to produce energy dependence 
in a damped field, as discussed in Section 2.5.2| 

Damping model fits to all other regions predict consis- 
tently thin rims with measurable FWHMs at decreasing 
energy. The width-energy dependence parameterized by 
Toe trends towards zero at low energy, indicating that 
rim widths are comparatively energy independent. At 
higher energy, such rims narrow slightly with increas- 
ing energy (ttze ~ —0.2) to match the observed width- 
energy dependence in X-rays. The gradual increase in 
|toe| with energy is expected in the damping model as 
advection and/or diffusion take control of rim widths at 
increasing energy, as discussed in Section [2~5l The best 
damped fits for these regions all have larger B 0 values 
than in the “weak-field” damping case. We refer to these 
fits as giving rise to “strong-field” damping, associated 
with stronger magnetic fields, weaker emission intensity 
behind the thin rim, and clear rims with measurable 
FWHMs at low photon energies (at and below soft X- 
rays). 

We are unable to determine whether weak- or strong- 
field damping is descriptive of Tycho’s shock magnetic 
field, given the large y 2 values on our model fits. Sev- 
eral regions may be well-fit by “weak-field” and “strong- 
field” damping alike. But, the qualitative behavior of Bq 
is suggestive; very strong fields at the shock (> 100 /zG) 
with magnetic field damping on a scale comparable to 
rim FWHMs appear incompatible with significant energy 
dependence. The distinction between weak- and strong- 
field damping may be equally well described as strong 
versus weak damping. The dichotomy reflects additional 
degeneracy between Of, and B 0 - for a given rim width, 
increased damping length ab requires increased Bq to re- 
produce the same width, and vice versa. 

Although model fits are poorly-constrained, we at- 
tempt to further explore how damping and synchrotron 
losses set rim widths for damped fits with finite dj,. 
Table [5] compares rim widths at 2 keV to advection 
lengthscale Z a d and damping lengthscale ab from the loss- 
limited and damped fits of Table [3j As discussed in 
Section |2.5[ loss-limited rim widths are qualitatively set 
by either damping or the dominant transport process as 
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Table 5 

Lengthscale analysis 


Region 

Measurements 

Loss-lim. fit 


Damped fit 


Vd 

(10 8 cm/s) 

w (2 keV) 
(%r 3 ) 

^ad 

(%r s ) 

Jad/^diff 

(-) 

^ad 

(%r s ) 

&b 

(%r s ) 

^ad/ 

(-) 

w ( 2 keV)/ab 
(-) 

1 

1.30 

2.64 

0.61 

1.40 

10.26 

0.80 

12.83 

3.30 

2 

1.29 

0.99 

0.27 

1.39 

11.54 

0.30 

38.48 

3.28 

3 

1.29 

0.74 

0.17 

1.38 

10.80 

0.20 

54.00 

3.72 

4 

1.28 

1.36 

0.31 

1.38 

10.08 

0.40 

25.20 

3.40 

5 

1.28 

1.27 

0.30 

1.37 

12.07 

0.30 

40.24 

4.25 

6 

1.28 

1.24 

0.18 

1.37 

2.61 

0.40 

6.53 

3.10 

7 

1.27 

0.96 

0.17 

1.36 

0.29 

0.60 

0.49 

1.61 

8 

1.27 

0.99 

0.19 

1.36 

0.90 

0.50 

1.80 

1.98 

9 

1.26 

1.03 

0.17 

1.36 

2.09 

0.40 

5.23 

2.57 

10 

1.26 

0.73 

0.14 

1.35 

9.21 

0.20 

46.05 

3.66 

11 

1.25 

1.05 

0.21 

1.34 

0.25 

1.00 

0.25 

1.05 

12 

1.24 

1.09 

0.25 

1.33 

10.80 

0.30 

35.98 

3.62 

13 

1.23 

0.98 

0.18 

1.32 

0.34 

0.60 

0.57 

1.64 

14 

1.14 

0.93 

0.17 

1.23 

2.84 

0.40 

7.10 

2.32 

15 

1.14 

0.75 

0.15 

1.22 

3.95 

0.30 

13.16 

2.50 

16 

1.13 

0.63 

0.12 

1.21 

0.14 

0.60 

0.24 

1.06 

17 

1.12 

0.64 

0.13 

1.20 

0.36 

0.40 

0.91 

1.61 

18 

1.15 

1.32 

0.28 

1.23 

11.29 

0.30 

37.64 

4.42 

19 

1.18 

0.95 

0.17 

1.27 

1.99 

0.40 

4.97 

2.37 

20 

1.17 

0.78 

0.14 

1.26 

3.08 

0.30 

10.28 

2.59 


Note. — All lengthscales computed at fiducial energy 2 keV from best fits with 77 2 = 1 and /x = 1. Ratio of / a d/^diff is same for loss-limited and 
damped fits and depends only on plasma velocity Vd and observation energy, as l a d/ldiff is independent of B 0 for fi = 1. 


w ~ min (ab, max (Z a( j, Zdiff))- At 2 keV with r ]2 = 1 and 
fi = 1 fixed, the ratio l ad /l dm is nearly constant; vari- 
ation (1.2-1. 4) arises solely from azimuthal variation in 
shock velocity. If rims are loss-limited and diffusion is 
negligible, we anticipate rim widths w = 4.6Z ad , where 
the factor 4.6 may be derived assuming spherical sym- 
metry and an exponential synchrotron emissivity (Ballet 


20061. At 2 keV, we find that loss-limited rim widths 


w ~ 5Z a <j due to diffusion. Enforcing w ( 2 keV)/db > 1 for 
damped fits still permits l ad /ab < 1 as w { 2 keV) < 5Z ad , 
with damping decreasing w ( 2 keV) below the expected 
loss-limited width at a given B 0 . We may impose a 
tighter or looser bound, but our results should not be 
greatly affected given large uncertainty in our fits and 
associated degeneracy between Bq and a&. 

Our modeling can also reassess the significance of rim 
width-energy dependence in the remnant of SN 1006 pre- 
sented by |R14[ We fit averaged filament widths in SN 
1006 measured by |R14| to our damping model using the 
same procedure as for Tycho. We fix B min = 5 p G, al- 
though best fit magnetic fields for SN 1006 are smaller 
than those of Tycho due to the much wider filaments 
of SN 1006. A more extensive search of parameter space 
produces acceptable damping (hybrid) models with more 
rapid shrinkage than found in R14[ Damped fits are com- 
parable to or better than loss- limited fits in 3 of 5 fila- 
ments in SN 1006. Fits to two filaments with strong en- 
ergy dependence (toe ~ —0.5) favor a loss-limited model 
with sub-Bohm diffusion ( rj 2 -C 1), though sub-Bohm 
diffusion may be an unphysical result pointing to over- 
simplifications in the model. The width-energy depen- 
dence in SN 1006 (toe ~ —0.3 to —0.5) is overall slightly 
stronger than in Tycho, and the best damped fits for SN 
1006 all fall into the “weak-field” case. The best fit Bq 
values are less than 40 / iG in the damped model, com- 
pared to ~ 100-200 fiG i n the loss-limited mode l. If thin 
radio rims in SN 1006 (Reynolds fc Gilmore||l986 ) indi- 


cate magnetic field damping on lengthscales comparable 
to X-ray filament widths, some magnetic damping is not 
incompatible with rim width narrowing. A future study 
of SN 1006 with recent multi-frequency Karl G. Jansky 
Very Large Array data (PI: D. Green) may verify whether 
radio and X-ray rims can be jointly described by an am- 
plified and subsequently damped magnetic field. 


5.2. Model assumptions and correctness 

Our models adopted a distance d to Tycho of 3 kpc, 
but es timates for Tycho’s distance range between 2.3- 
4 kpc (|Hayato et al. 2010). As a larger remnant distance 
would increase both physical filament widths and shock 
velocity estimates from proper motion, we may derive 
expected scalings for modeled 772 and B 0 as functions of 
assumed distance d. The advective lengthscale, equa- 


tion (10 1, may be rearranged to obtain: 


B 0 = (3.17 pG) 
lad 


Vd 

10 s cm/s 
-2/3 


0.01 kpc 


2/3 


hv 

1 keV 


-1/3 


or, more simply, 

Bq oc ( Vd ) 21 3 (? ad ) _2/3 W 1/3 . 


(15) 


(16) 


Both l a d and Vd scale linearly with remnant distance d 
and thus their effects cancel in determining the mag- 
netic field. If diffusion is th e primary control on filament 
lengthscales, equation (fTTl) yields: 


V2 0c(l diS ) 2 B^ +5)/2 v-^~ 1 V 2 


(17) 


Model fits with varying distance obey both scalings, 
r \ 2 oc d 2 and Bq constant. When comparing model fits 
with remnant distances of 3 kpc and 4 kpc, the devia- 
tion from the idealized scaling is < 1% for Bq and ~1- 
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5% for 772 - As varying remnant distance d leaves width- 
energy scaling ttie invariant, the relative contributions of 
/ a d and Zdiff should also be invariant. Then both length- 
scales should scale simultaneously with d , yielding the 
observed behavior. In the damped model, a larger dis- 
tance d will require larger physical damping lengths from 
magnetic turbulence. But, fitted values of ab should re- 
main unchanged as we report ab in units of shock radius 

r s - 

The exponential variation in D[x) (equation (16])) may 
not be physically reasonable; as noted above, the as- 
sumption D(x) oc 1 /B 2 [x) gives rise to sharp gradients 
in diffusion coefficient and even inverts the effect of 77 2 on 
our model profiles (where larger 772 can cause rim widths 
to narrow). Nevertheless, the modeled behavior appears 
physically reasonable. Only X-ray profiles are impacted 
by the assumption on D(x ) as radio profiles assume no 
diffusion. And, model behavior driven by advection and 
magnetic damping - namely, rim width-energy depen- 
dence in a damped magnetic field - occurs for D = 0, and 
cannot be an artifact of our assumption that D{x)B 2 (x) 
is constant. 

We emphasize that additional counts from averaging 
measurements or selecting larger regions will likely not 
improve our ability to constrain Bq and 772 from width- 
energy modeling. This is easily seen from Figure[5]as well 
as large y 2 values in Tables [ 3 ] that reflect relatively tight 
errors in our FWHM profile measurements (Table [ 2 ]). 


5.3. Other potential constraints 

Rim width-energy dependence is a morphological man- 
ifestation of spectral softening downstream of the for- 
ward shock; previously, Cassam-Chenai et al. (2007) also 
sought to distinguish loss-limited and damped rims with 
a careful spectral study and 1-D hydrodynamical model. 
Although knowledge of radial spectral variation, in prin- 
ciple, fully determines rim profiles, there is not a clear 
relationship between observed spectral variation and rim 
widths alone. Our model, for example, replicates ob- 
served rim width variation but underpredicts the ob- 
served s pectral variation . The more sophisticated model 
of Cassam-Chenai et al. (2007) similarly has trouble re- 
producing the observed radial gradient in spectral pho- 
ton index. Photon indices predicted in our model, for 
integrated model rim spectra similar to those of Figure [3] 
and Table [T] are also somewhat ill-constrained. Enforc- 
ing limi ts on acceptable model photon indices (e.g., Fig. 
16 Cassam-Chenai et al. 2007) could help constrain shock 
parameters or identity discrepancies between model as- 
sumptions and measurements. 

Estimates of magnetic field strength that depend upon 


global models for Tycho’s evolution (e.g., Morlino & 
Caprioli 2012 ) may not help constrain our rim model 


results. In such models, assumptions on particle accel- 
eration and magnetic field evolution downstream of the 
shock (e.g., further damping or amplification at the re- 
verse shock and contact discontinuity) would affect the 
spatially integrated spectrum of Tycho. We have focused 
on effects immediately behind the forward shock, making 
as few assumptions as possible about particle accelera- 
tion and magnetic fields throughout the remnant. It may 
not make sense to extrapolate our results to estimate 
magnetic fields throughout the remnant or, similarly, to 
use global models to constrain our results given all as- 



Figure 7. Radio image of Tycho’s SNR at 1.375 GHz with linear 
scaling. Extraction regions (green) for joint radio and X-ray pro- 
file analysis overlay region selections for X-ray rim width analysis 
(Figure IT] . 


sumptions involved. 


6. JOINT RADIO AND X-RAY MODELING 

Thin radio rims spatially coincident with X-ray syn- 
chrotron rims may help constrain magnetic field damp- 
ing. Radio synchrotron emission in the remnant interior 
may arise from not only recently shock-accelerated elec- 
trons but also long-lived electrons interacting with a tur- 
bulent field inside the remnant, as ~GeV electrons have 
cooling times of order 10 5 -10 7 yr in 10-100 /xG magnetic 
fields. But simple models for synchrotron emissivity as 
a function of density in Sedov and pre-Sedov dynami- 
cal stages predict radio synchrotron profiles rising grad- 
ually to broad maxima well inside the shoc k radius (e.g. 
Reynolds & Chevalier|1981||Reynolds|1988|). Our steady- 
state planar transport model predicts monotonically in- 
creasing emission downstream of the forward shock in 
a loss-limited model, until sphericity becomes important 
and our model is inapplicable. On scales of a few percent 
of the remnant radius, only magnetic field damping can 
cause emission to decrease. 

The idea of modelin g X-ray and radio pr ofiles jointly 
was pioneered by |Cassam-Chena'f et al.| (|2007|), who 
found that sharp radio rims in lycKo were not repro- 
duced by a loss-limited model. Cassam-Chenai et al.| 
(2007) considered loss-limited and damping profiles con- 


sistent with physical constraints applied in a hydrody- 
namical model. Here, we neglect physical constraints 
and estimate a damping length necessary to generate ra- 
dio rims with shape similar to observations. 

We extract radio profiles from a 1.375 GHz image of 
Tycho taken with the Very Large Array (VLA) in A con- 


figuration in March 1994 (PI: D. Moffett); see Reynoso 
et al. (1997) for a detailed presentation. The halt-power 
beam width of ~1.5" just resolves thin radio rims and 
structure near the forward shock; the image is sampled 
at 0.5". We also extracted 4-7 keV X-ray profiles in 
all regions from the previous archival Chandra observa- 
tion to jointly model radio and X-ray profiles, permit- 
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Table 6 

Fit parameters from profile shape comparison 


Region 

B 0 

(M g) 

ab 

(-) 

Region 

B 0 

(MG) 

Ub 

(-) 

A 

50 

0.020 

I 

250 

0.020 

B 

200 

0.050 

j 

300 

oo 

c 

15 

0.010 

K 

400 

0.010 

D 

100 

0.050 

L 

250 

oo 

E 

120 

0.030 

M 

200 

oo 

F 

300 

0.025 

N 

800 

0.010 

G 

250 

0.020 

o 

200 

0.005 

H 

300 

0.020 

p 

150 

0.012 


Note. — Damping lengths of oo indicate that a loss-limited fit is 
favored (a& > 10% of shock radius r s ). 


ting somewhat firmer discrimination of plausible model 
parameters. Figure [7] shows the extraction regions over- 
laying the radio image and the previous X-ray profile 
regions of Figure [I] 

We compute model radio and X-ray measured profiles 
from the transport model of equation |l]) for varying Bq 
and a,b, similar to those shown in Figure [2] The pa- 
rameters B m i n = 5 yuG, 772 = 1, and p. = 1 are held 
fixed. Diffusion, in particular, is negligible for modeled 
radio emission as particle energies are 3 orders of mag- 
nitude lower than in X-ray. Neglecting diffusion also cir- 
cumvents the unphysical assumption that D(x)B 2 (x) is 
constant, which was invoked to obtain Green’s function 
solutions to equation ([I]). 

We align each set of model profiles at some B 0 and 
db to measurements by eye, varying relative amplitudes 
and translations in radio and X-ray independently to 
best match the measured profiles. Although “fitting” 
by eye does not quantitatively bound model parame- 
ters, we can find plausible values of Bq and a& and es- 
timate the importance of magnetic damping throughout 
the remnant. More involved nonlinear fitting may be un- 
reliable as we cannot constrain spatially heterogeneous 
radio emission within the remnant. Moreover, our trans- 
port model neglects self-similar downstream evolution of 
shocked plasma (decaying velocity, density, and magnetic 
field) and is inaccurate further downstream than ~10% 
of the shock radius. 

The joint radio and X-ray profile modeling contrasts 
strongly with our previous width-energy fitting from X- 
ray measurements alone, which neglected profile shape in 
favor of more robust FWHM measurements. Manually 
fitting profiles allows us to consider radio and X-ray fila- 
ments that do not have well-defined FWHMs, especially 
as radio rims do not fall below 50% of the peak emission. 
We can also use profiles from regions not previously con- 
sidered due to the lack of an X-ray FWHM, especially in 
softer (0.7-4 keV) X-rays. 

We identify three classes of radio profiles: thin rims 
with downstream troughs, plateaus, and continuous rises. 
Regions B, C, and D (southern limb) show plateaus in 
radio emission. Regions J, L, and M (around NW) show 
continuous rises in emission. All other regions have a 
radio rim within 15" of the forward shock, where the for- 
ward shock in radio is assumed at zero intensity. Figure [8] 
shows extracted radio and nonthermal X-ray (4-7 keVj 
profiles for three well-fit regions with our best manually 
selected model profiles, illustrating each of the three ra- 
dio profile types observed. 
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Figure 8. Measured radio and X-ray profiles plotted with model 
profiles for varying a 5 and Bo in each region, showing typical pa- 
rameters (and ranges) required to reproduce radio and X-ray rims 
simultaneously in our model. Solid black curves in all regions plot 
our manually chosen best model profiles. Profiles are chosen to 
show varying radio rim morphology, including plateaus (B), thin 
rims with troughs (K), and continuous rises (M); these profiles 
show some of the best agreement of our selected regions, but cf. 
Figure [9] Profile radial coordinates are shifted arbitrarily to aid 
visual comparison. Negative radio intensity is unphysical and as- 
sociated with deconvolution of raw VLA visibilities. 


Our model requires damping length ab < 0.1 to pro- 
duce a plateau or thin rim in radio emission. For regions 
with thin radio rims, the best manually selected profiles 
have B 0 between 50-400 fi G, neglecting only Region N 
which could not be modeled simultaneously in both X- 
ray and radio (Figure^). The damping length ab ranges 
between 0.01-0.03, or2 7" . We list estimated best fit 
parameters in Table [6] 

The radio plateaus in regions B, C, D are compatible 
with damping lengths between 0.01-0.05. Continuous 
rises in radio emission are best modeled with a& > 0.1 and 
l?o~200-300 ^iG. Although damping lengths a& > 0.1 
are not physically meaningful well beyond the shock, 
these large values of ab yield practically constant mag- 
netic field near the shock. Our results are also insen- 
sitive to the assumed value of R m i n ; only model pro- 
files with small magnetic fields (Bq < 100 /rG) and/or 
small damping length (ab < 0.01) are affected if we take 
B m in = 2 /i G rather than 5 /rG. 

Tycho’s shock structure is more complex than assumed 
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Figure 9. Extracted profiles poorly reproduced by our model, 
compared to Figure [8] Region I shows irregularly shaped radio 
rim. Region N contains two superposed filaments that cannot be 
modeled by a single rim in both X-ray and radio; the narrow X- 
ray rim requires atypically strong magnetic fields (~ 800 /iG), and 
the emission plateau behind the radio rim requires small damping 
lengths (a& ~ 0.005). 


by our transport model. Emission towards the remnant 
interior clearly shows spatial structure (Figure[7|. Figure 
[9] shows two regions that were poorly described by our 
model. A majority of our regions have irregular rims; 
in at least 2-3 regions, this may be attributed to pro- 
jection of multiple filaments. Others (e.g., Region K) 
show rims with slopes that cannot be matched by our 
models, whether too steep or shallow. Shape mismatch 
may be attributed in part to point-spread mismatch, dif- 
fusion (e.g., 772 7 ^ 1 ), shock precursors, or other effects. 
Nevertheless, the conclusion that thin radio rims require 
magnetic damping is supported by more sophisticated 
modeling. As we have noted, h ydro d ynamic mo dels with 
diffu sive shock acceleration (Cassam-Chenai et al.| 2007 


Slane et al.|2014 ) also cannot produce radio profiles with 


narrow rims in a purely advected magnetic field. It is 
also not clear that radio structure is due to a globa l ra- 
di al va riation in magnetic field strength. E.g., [Slane et ah] 
(2014) suggest that Rayleigh- Taylor fingers between the 
forward shock and contact discontinuity could locally 
confine radio-emitting electrons, creating observed radio 
rim structure. But, if radio rims are due strictly to mag- 
netic field effects, our conclusions on magnetic field drop- 
off are independent of the causative physical mechanism 
so long as they yield magnetic field fall-off over 1 % of 
remnant radius. 

We may also use spatially overlapping radio and X- 
ray region selections (Figure [7]) to attempt to constrain 
X-rav width-energy fits. All OUray region selections (Fig- 
ured]) are associated with radio rims, except for Regions 
11 and 12 underlying Region J in radio. X-ray emission 
at regions B, C, D (radio plateaus) and L, M (radio rises) 
showed rims, but their widths either could not be mea- 
sured, or could only be measured in the 4-7 keV band. 


If radio rims require damping at ~1% remnant radius, 
width-energy fits with “weak-field” damping lengthscales 
of ~0.5% may be disfavored. Conversely, the radio and 
X-ray filament of Regions 11,12 and J is the best remain- 
ing loss-limited rim candidate, although width-energy fits 
are equivocal towards damped and loss-limited models. 


7. CONCLUSIONS 

We measured the widths of several thin synchrotron 
filaments around Tycho’s supernova remnant and found 
moderate narrowing of rim widths throughout the rem- 
n ant, corro borating rim narrowing observed by |Ressler| 
et al. (2014) in the remnant of SN 1006. We confirmed 
that selected filaments are dominated by nonthermal 
emission and have clearly measurable full widths at half 
maximum in 4-5 energy bands. Both X-ray width-energy 
fits and joint radio/X-ray profile modeling require mag- 
netic fields > 20 /rG even with magnetic damping. 

A steady-state particle transport model with constant 
magnetic field gives diffusion coefficients and magnetic 
field strengths broadly consistent with prior estimates 
from rim widths (e.g., Parizot et al. 2006 Rettig fe 


Pohl|20^|) and r adio ana gamma ray m easurements ( Ac- 


ciari et all 2011 |Morlino fe Caprioli [2012 ) . The same 
model with a damped magnetic held is equally capa- 
ble of describing our measured data. At weak energy 
dependence the two models are indistinguishable and 
magnetic damping fits favor moderately amplified mag- 
netic fields beyond simple compression, but lower than 
for loss-limited models. At moderate energy dependence 
(toe ~ —0.3), the damping model permits weak magnetic 
fields and short damping lengths (< 1 % remnant radius) 
to reproduce energy dependence, but we still cannot fa- 
vor either damped or loss-limited rims due to fit uncer- 
tainty. The distinction between loss-limited and damped 
models is somewhat artificial; only for large differences 
between the various transport lengthscales, a frequency- 
dependent occurrence, is one or the other mechanism 
clearly dominant. 

Thin radio synchrotron rims, however, are not re- 


produ ced in loss-limited models (Cassam-Chenai et al. 
2007 ) . Assuming shocked electrons account for most ra- 


dio emission immediately downstream of Tycho’s forward 
shock, we jointly model radio and X-ray profiles and find 
that damping lengths of 1 -5% of the shock radius are 
required throughout most of the remnant; only a few 
(3/16) selected regions are plausibly consistent with a 
constant advected magnetic field. Typical magnetic field 
strengths range between 50-400 /iG. Although we can- 
not bound damping lengths and fields from our qualita- 
tive profile comparisons, our results are physically rea- 
sonable and are likely good to order-of-magnitude. If 
damping lengths inferred from radio rims are correct, 
“weak-field” damping is disfavored in explaining X-ray 
rim width-energy dependence, and damped rims require 
magnetic field amplification to ~100 7 iG or more in Ty- 
cho. 


We thank the anonymous referee for comments that 
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